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Density- functional theory calculations based on conventional as well as hybrid exchange-correlation functionals 
have been carried out to study the properties of helium in various oxides (AI2O3, Ti02, Y2O3, YAP, YAG, 
YAM, MgO, CaO, BaO, SrO) as well as at oxide-iron interfaces. Helium interstitials in bulk oxides are 
shown to be energetically more favorable than substitutional helium, yet helium binds to existing vacancies. 
The solubility of He in oxides is systematically higher than in iron and scales with the free volume at the 
interstitial site nearly independently of the chemical composition of the oxide. In most oxides He migration 
is significantly slower and He-He binding is much weaker than in iron. To quantify the solubility of helium 
at oxide-iron interfaces two prototypical systems are considered (FejMgO, Fe|FeO|MgO). In both cases the 
He solubility is markedly enhanced in the interface compared to either of the bulk phases. The results of the 
calculations allow to construct a schematic energy landscape for He interstitials in iron. The implications 
of these results are discussed in the context of helium sequestration in oxide dispersion strengthened steels, 
including the effects of interfaces and lattice strain. 
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I. INTRODUCTION 

Fusion environments are characterized by an abun- 
dance of energetic helium ions (a-particles, He^+) that 
are produced in the fusion reaction as well as via 
nuclear transmutation reactions in the plasma facing 
components. '^'^ Candidate materials for structural appli- 
cations in fusion reactors must therefore be able to tol- 
erate large concentrations of helium without degradation 
of their mechanical properties under extreme conditions 
of up to 200 displacements per atom and 2000 appm He.^ 
Most conventional steels are unable to sustain these con- 
ditions since they suffer from He mediated bubble as well 
as void formation and growth ( "swelling" ) leading to em- 
brittlement and mechanical failure.'^ 

Using first-principles calculations it has been shown 
that in iron (as a model system for steels) He prefers 
to occupy substitutional rather than interstitial sites. ^'^ 
For dynamical reasons, however, during irradiation the 
majority of He is introduced in the form of interstitials.^ 
Since the effective interaction between He interstitials in 
the iron matrix is attractive^ they can readily form clus- 
ters. The latter are further stabilized by the addition of 
vacancies, which eventually leads to the formation and 
growth of bubbles and voids. 

To overcome the shortcomings of conventional steels 
for applications in fusion environments, it has been 
suggested to employ nano-structured ferritic alloys,'^ in 
particular oxide dispersion strengthened (ODS) steels. 
These materials are obtained by mechanical alloying of 
oxide particles with steel powders (see e.g., Refs. 7- 
9). They are characterized by a fine distribution of 
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nanometer-sized oxide particles"'^^""'^^ that act as obsta- 
cles for dislocation motion and are metastable up to very 
high temperatures. The oxide particles and even more so 
the oxide-matrix interfaces are expected to act as sinks 
for He interstitials. The very high density of these par- 
ticles should lead to a fine distribution of He bubbles 
and thereby effectively limit the formation of larger su- 
percritical voids that lead to mechanical failure.'^ Sev- 
eral experiments provided evidence that ODS steels are 
much more swelling resistant than their non-oxide con- 
taining counterparts, raising hopes that these materials 
will eventually satisfy the design criteria for future fusion 
reactors. 

Designing ODS steels that can sustain the extreme con- 
ditions inside a fusion reactor requires a close collabora- 
tion between experiment and modeling. As there are 
currently no neutron sources available that can repro- 
duce the intense neutron spectrum resulting from fusion, 
a combination of (experimental and numerical) simula- 
tions must be employed to predict material performance. 
In this context, numerical modeling of the time evolu- 
tion of defect populations (vacancies, interstitials, inter- 
stitial loops, dislocations, bubbles, voids etc.) plays a 
pivotal role.-'^'^'-'^^ Such models rely on a database of rates 
for various microscopic processes that occur in the ma- 
terial, which typically comprises both experimental and 
atomic scale modeling data. Given the importance of 
oxides particles in improving swelling resistance, the mi- 
croscopic mechanisms that govern their interaction with 
He deserve particular attention, yet at present our mi- 
croscopic understanding of these interactions is limited. 

A recent very extensive transmission electron mi- 
croscopy study^^ revealed that in ODS steels bubbles 
and voids form preferentially in the close vicinity of ox- 
ide particles, leading in many cases to the formation of a 
"ring" of bubbles surrounding a particle. The same study 
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also demonstrated that oxide particles in ODS steels ex- 
hibit a broad variety of size, chemical composition, bulk 
and interface structure including amorphous, partially 
amorphous and crystalline particles, core-shell configu- 
rations as well as chemical gradients. It was further- 
more shown that materials with smaller particles can 
tolerate more helium providing direct evidence for the 
importance of the interface area.^^ In this context, an 
oxide "particle" can be as small as a few Angstroms, 
correponding to just a couple of atoms. While size of 
these "nanofeatures" ^ can be below the resolution obtain- 
able in transmission electron microscopy, they nonethe- 
less contribute to He sequestration due to their sheer 
number and very large effcetive interface area. 

Atomic scale modeling is in principle ideally suited 
to complement these experiments and to provide de- 
tailed insight as well as quantitative information regard- 
ing the behavior of helium inside and near oxide parti- 
cles. The enormous chemical and structural complexity 
present in ODS steels, however, renders a direct simu- 
lation of these systems at present impractical. On the 
other hand macroscopic measurements of He populations 
indicate that the most important parameter is particle; 
size and thus interface area.^^ This suggests that the sol- 
ubility of helium follows more general trends that are 
independent of variations in the local chemistry. The ob- 
jective of the present paper is to demonstrate that the 
behavior of He in oxides and in the vicinity of oxide-iron 
interfaces can to a large extent be described by simple 
scaling relations. This is achieved by means of density- 
functional theory calculations using both conventional 
and hybrid exchange-correlation (XC) functional of bulk 
oxides, bulk iron, and representative oxide-iron inter- 
faces. The thus obtained data provides not only valu- 
able insight into the microscopic mechanisms but can be 
further utilized to parametrize for example rate equation 
models.^'* Note that even though nanosized oxide nuclei 
of the type described in Refs. 10 and 12 are not explicitly 
studied in the present work, the arguments that indicate 
that the governing parameter is free volume and in par- 
ticular free volume at oxide- iron interfaces in general also 
transpire to the case of very small oxide inclusion. Com- 
pared to earlier investigations that considered He defects 
in select oxides and carbides, the present work aims 
to provide a more general perspective, including a vari- 
ety of oxides that represent different local environments, 
chemistry, and covalent-ws- ionic character. 

The paper is organized as follows. In the next section 
computational methods and parameters are summarized. 
The results of a comprehensive study of He-related de- 
fects in three prototypical oxides (AI2O3, Ti02, Y2O3) 
are described in Sect. HI. It is demonstrated that He oc- 
curs predominantly in the form of interstitials and that 
solubilities are much larger than in iron. Furthermore 
it is found that formation energies of He interstitials 
determined using conventional XC functionals are very 
close to values obtained from (more elaborate) hybrid XC 
functional calculations. To rationalize the variation of 



He interstitial formation energies among different oxides, 
in Sect. IV A their volume dependence is studied addi- 
tionally including the oxides Y4AI2O9 (YAM), Y3AI5O12 
(YAG), YAIO3 (YAP), MgO, CaO, SrO, and BaO. Sec- 
tion IV B concerns the migration barriers for He inter- 
stitials in oxides, which are found to be systematically 
higher than in iron. The binding between He interstitials 
is the subject of Sect. IV C, where it is shown that most 
oxides are less prone to He cluster formation than Fe, 
a behavior that results from a lower density of intersti- 
tial site in these materials. The solubility of He intersti- 
tials in two representative oxide- iron interfaces (FejMgO, 
Fc|FeO|MgO) is quantified in Sect. V, where solubilities 
at interfaces are found to be systematically higher than 
in the bulk phases. Finally all these data are combined 
to sketch a typical energy landscape for He interstitial 
migration across an oxide-iron interface and discuss the 
results in the context of He sequestration in ODS steels. 



II. METHODOLOGY 

A. Thermodynamics 

Whereas in elemental metals the formation energy of 
an intrinsic defect is constant, in oxides defect formation 
energies depend both on the electron chemical potential 
(also referred to as Fermi level) and the chemical envi- 
ronment. The expression for the defect formation energy 
makes these dependencies explicit, 

AEd = {Ed - Eh) - qiEyBM + Me) - X] ^^H^i' (1) 

i 

where /i^ is the chemical potential of component i, Ejj 
is the total energy of the system containing the defect 
and Eh is the total energy of the ideal reference system. 
The defect formation energy is linearly dependent on the 
defect charge state q and the electron chemical potential 
/Ue which is measured with respect to the valence band 
maximum EysM- Here, Arij denotes the difference in 
the number of atoms of element i between the system 
with and without the defect, for example in the case of 
an isolated oxygen vacancy Ano = — 1 whereas all the 
other Arii are zero. The chemical potentials of the con- 
stituents /i; are conveniently expressed with respect to 
the bulk chemical potentials /i^""^ of the respective el- 
emental ground states, /Xj = /x^"'"^ + A^ii. Metal and 
oxygen-rich conditions correspond to A/imctai = and 
A^o = 0, respectively. The values of the chemical po- 
tentials can be translated to partial pressures enabling 
direct comparison with experiments.^^ For simplicity, in 
the present work He-rich conditions are assumed always, 
i.e. A/iHc = cV. For an elemental metal Eq. (1) reduces 
to the usual expression 

AEd = Ed r^^ — Eh, (2) 
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TABLE I. Overview of computational parameters used in calculations of properties of the ideal bulk systems as well as point 
defects. Migration barrier calculations for Y2O3 and the rocksalt structured oxides were carried using the parameters given in 
brackets. 



A1203 


Ti02 


Y2O3 


YAP 


YAG 


YAM 


(Mg,Ca,Ba,Sr)0 


Ideal cell calculations 














Number of atoms 10 


6 


40 


20 


80 


60 


2 


fc-point sampling r4 x 4 x 4 


re X 6 X 6 


r2 X 2 X 2 


r4 X 4 X 4 


r2 X 2 X 2 


r2 X 2 X 2 


FS X 8 X 8 


Defect calculations 














Number of atoms 270 


216 


320" 


160 


160 


60 


216** 


Type of supercell rhombohedral 


tetragonal 


body-centered cubic orthorhombic simple cubic monoclinic 


simple cubic 


Supercell size 3x3x3 


3x3x4 


2x2x2 


2x2x2 


1x1x1 


1x1x1 


3x3x3 


fc-point sampling Fl x 1 x 1 


n X 1 X 1 


n X 1 X 1 


ri X 1 X 1 


ri X 1 X 1 


ri X 1 X 1 


ri X 1 X 1 


XC functional PBE, HSE06 


LDA, HSE06 


PBE 


PBE 


PBE 


PBE 


PBE 



" Migration barriers calculated using 80-atom cells and a r2 x 2 x 2 fe-point sampling. 
Migration barriers calculated using 64-atom cells and a r4 x 4 x 4 fc-point sampling. 



where A'^ denotes the number of atoms in the ideal cell. 

In terms of their formation energies, the binding energy 
of two defects A and B is given by 

i\E^{AB) = AEf{AB) - AEf{A) - AEf{B). (3) 

Following this convention negative binding energies cor- 
respond to exothermic defect reactions and imply an at- 
tractive interaction between A and B. Note that the 
binding energy can change as a function of the chemical 
potential but is independent of the chemical environment 
since the An^ terms that appear in Eq. (1) cancel each 
other in Eq. (3). 

Assuming independent defects (low density limit the 
equilibrium concentration of a defect is related to its free 
energy of formation AGf according to 

Ceg = Co exp [-AGf/kBT] , (4) 

where the pre-factor cq denotes the density of potential 
defect sites per unit volume. For extrinsic defects the 
equilibrium concentration corresponds to the solubility 
of the foreign element in the matrix. The free energy of 
defect formation AGf can be decomposed into the for- 
mation energy AEf as well as the vibrational TASf^mb 
and electronic TASf^ei formation entropies. (Note that 
Eq. (4) already incorporates the configurational entropy, 
see Ref. 22). For materials with a band gap {Eg » ksT) 
the electronic contribution to the formation entropy is 
virtually zero and even for metals this term is usually 
small compared to the other contributions. The vibra- 
tional formation entropy can become important at el- 
evated temperatures. The present work is, however, 
mainly concerned with comparing defects with very sim- 
ilar characteristics and therefore relative changes of the 
formation entropy between different defects can be ex- 
pected to be small compared to the formation energy 
term. In the following, therefore only formation energies 
are considered and AGf ~ AEf. 

It should be stressed that formation energies are equi- 
librium quantities while a material under intense irradi- 
ation is obviously a non- equilibrium system. Yet within 



the constraints of such a scenario, it is nonetheless useful 
to consider formation energies as they will determine the 
driving forces, which determine the long time evolution 
of the system. 



B. Computational details 

Calculations were carried out within density- 
functional theory using the projector augmented 
wave formalism^^'^^ as implemented in the Vienna 
ab-initio simulation package. Ti-3p, Y-4s, Y-4p, 
Ca-3p, Sr-4s, Sr-4p, Ba-5s as well as Ba-5p states were 
treated as part of the valence. The plane-wave energy 
cutoff was set to 500 eV for all calculations. To represent 
exchange and correlation effects, we employed the local 
density approximation (LDA), the generalized gradient 
approximation as parametrized by Perdew, Burke and 
Ernzerhof (PBE)^^ as well as range-separated hybrid 
functionals^*^ obtained by mixing conventional XC 
functionals (LDA or PBE) with 25% exact exchange 
at short-range with a screening parameter of 0.2 A . 
These functionals are refered to as HSE06lda and 
HSE06, respectively. 

For LDA and PBE the electronic contributions to the 
dielectric constants reported in Table III were computed 
within the linear response approach taking into account 
local field effects. For the hybrid functionals the elec- 
tronic contribution was calculated from matrix elements 
of the dipolc operator in the velocity gauge and the lo- 
cal field effect correction from either LDA or PBE was 
added. The ionic contribution was computed for LDA 
and PBE using linear response theory. 

Details regarding the defect calculations in oxides, 
specifically Brillouin zone sampling, supercell cell shapes 
and sizes, are summarized in Table I. For pure iron the 
PBE functional was used as well as 4 x 4 x 4 supercells 
containing 128 atoms and a 3 x 3 x 3 Monkhorst-Pack 
grid for sampling the Brillouin zone. For charged defects 
the monopole-monopole correction according to Makov 
and Payne'^^ was applied using the calculated static di- 
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electric constants given in Table III. Migration barri- 
ers were obtained via the climbing image-nudged clas- 
tic band method^^'^^ using three intermediate images to 
represent the transition path. The convergence of the 
calculations with respect to Brillouin zone sampling and 
supercell size was carefully tested, based on which the 
error in the He interstitial formation energies due to the 
computational parameters is estimated to be less than 
0.1 cV. 

Interfaces were modeled using slab geometries employ- 
ing a similar approach as in Ref. 34. For the ideal Fe|MgO 
interface the supercell contained six Fe and six MgO lay- 
ers equivalent to 36 atoms. The Fe|FeO|MgO interface 
model was composed of four Fe, two FeO, and six MgO 
layers equivalent to 40 atoms. Both systems were fully 
relaxed including cell shape and volume until forces were 
below 20meV/A and the components of the stress ten- 
sor less than Ikbar. For calculations involving He in- 
terstitials the supercell size was doubled parallel to the 
interface leading to supercclls with 144 and 160 atoms, 
respectively. Helium interstitial positions were system- 
atically sampled across the interface as well as the two 
"bulk" parts corresponding to 18 distinct configurations 
for each interface model. Each defect configuration was 
relaxed at fixed cell shape and volume until the maximum 
force fell below 30meV/A. In these calculations the Bril- 
louin zone was sampled using alOxlOxl Monkhorst- 
Pack grid for the ideal cells and a 5 x 5 x 1 grid for the 
defect cells. 

Results for bulk oxides considered in this study and 
an extensive comparison of different XC functionals is 
provided in the appendix. 

III. HELIUM AND INTRINSIC DEFECTS 

Helium can be incorporated either substitutionally or 
as an interstitial defect. In the next three sections, the 
thermodynamics of these two forms of He in three dif- 
ferent prototypical oxides (AI2O3, Ti02, Y2O3) arc com- 
pared. It is shown that interstitial He is the most stable 
form under most conditions and is also the most relevant 
form with regard to sequestration in ODS steels. 

A. Alumina 

Figure 1 shows the formation energies of interstitial 
and substitutional He as well as several intrinsic point de- 
fects in alumina. According to the calculations He inter- 
stitials preferentially occupy positions that are equivalent 
to 2b Wyckoff sites of the ideal structure. The interstitial 
formation energy is independent of the electron chemical 
potential, which is expected based on the closed-shell na- 
ture of He, and affected only slightly by the choice of XC 
functional (PBE: 2.15eV, HSE06: 2.23eV). 

When He is substituted for Al the formation energy of 
the resulting defect closely follows the Fermi level depen- 



dence of the Al vacancy safe for a constant upward shift. 
In contrast, while substitutional He on an oxygen site 
exhibits the same charge states as the oxygen vacancy, 
the 2-t- /O equilibrium transition level for Heo is strongly 
shifted with respect to Vq. 

This behavior can be rationalized by considering the 
relaxation patterns of different defect charge states. 
Since He has a closed shell the formation energy differ- 
ence between a vacancy and the corresponding substi- 
tutional defect results predominantly from strain. For 
the aluminum vacancy all charge states exhibit the same 
type of outward relaxation pattern and the formation 
energy difference between Vai and Hcai is only weakly 
affected by the charge state. In the case of the oxygen 
vacancy, however, the first neighbor shell relaxes inward 
for the neutral and outward for the positive (2+) charge 
state. This behavior is typical for "deep" oxygen vacan- 
cies and also observed in other oxides including e.g., ZnO 
and In203 (see Ref. 35).^^ As a result the strain energy 
contribution to the formation energies of Heo is strongly 
charge state dependent, which explains the strong shift 
of the equilibrium transition levels from Vq to Heo that 
is observable in Fig. 1. 

The absolute values of vacancy and substitutional He 
formation energies depend on the treatment of XC. The 
HSE06 hybrid functional, which provides a much im- 
proved value for the band gap compared to the PBE func- 
tional (compare Table HI), yields larger formation ener- 
gies approximately in accord with the increase in band 
gap. Yet PBE and HSE06 yield qualitatively the same 
picture with the equilibrium transition levels of Vq /Hcq 
and Vai/Hcai tracking the conduction and valence band 
edges, respectively. 

Figures l(b,c) illustrate not only the dependence of de- 
fect formation energies on the electron chemical potential 
but also on the chemical environment (compare Eq. (1) 
and discussion thereafter). Moving from Al-rich to O- 
rich conditions, which can be achieved by regulating the 
oxygen partial pressure, shifts the balance between Heo 
(Ano = — 1) and HeAi (AriAi = — 1) but leaves the for- 
mation energy of He^ unchanged (AriAi = Ano = 0). 
Figure 1 shows that the formation energies of both va- 
cancies and substitutional He can become negative under 
certain conditions. As negative formation energies imply 
the material being unstable with respect to defect for- 
mation, the Fermi level is constrained to the range in 
which all formation energies of intrinsic defects are posi- 
tive. Considering the thermodynamically allowed regions 
in Fig. 1, it can be concluded that for most conditions He 
interstitials are thermodynamically the preferred form of 
He in alumina. Yet He interstitials do bind to existing 
vacancies, since the reaction 

He, + Vx Hex 

is exothermic and the binding energy Eb^Rex) [compare 
Eq. (3)] is negative as shown in Fig. l(d,e). The figure 
also demonstrates that the impact of the XC functional 
on the binding energies is small. 
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Electron chemical potential (eV) Electron chemical potential (eV) Electron chemical potential (eV) Electron chemical potential (eV) 



FIG. 1. (a-c) Formation energies for point defects in alumina (a) calculated using the PBE XC functional under Al-rich 
conditions, and calculated using the HSE06 hybrid functional under (b) Al-rich and (c) O-rich conditions, respectively. The 
line slopes correspond to the defect charge states according to Eq. (1). (d-e) Binding energies of He interstitials to vacancies 
for (d) oxygen and (e) aluminum. 




FIG. 2. Formation energies for point defects in titania calcu- 
lated using (a) the PBE XC functional under Al-rich condi- 
tions as well as the HSE06 hybrid functional under (b) Ti-rich 
and (c) O-rich conditions. 



FIG. 3. Formation energies for point defects in yttria calcu- 
lated using the PBE XC functional under (a) Y-rich and (b) 
O-rich conditions, (c) Binding energies of He interstitials to 
vacancies. 



B. Titania 

The formation energies of several forms of He in tita- 
nia as well as related intrinsic defects are shown in Fig. 2. 
The lowest interstitial formation energy is obtained if He 
occupies positions equivalent to the 4c Wyckoff site of 
the ideal structure. Similar to the case of alumina sub- 
stitution on the cation sublattice leads to a defect with 
characteristics that are very similar to the cation vacancy. 
Unlike alumina this resemblance is also observed for oxy- 
gen vacancy and Heo. Once again this behavior can be 
related to the relaxation patterns of the oxygen vacancy. 
In contrast to alumina, the oxygen vacancy in titania is 
"shallow" , relaxation occurs inward for all charge states, 
and the strain energy associated with He insertion at 
the vacant oxygen site is virtually independent of charge 
state. The binding energy between a He interstitial and 
a vacancy is -l.OeV (LDA)/-l.leV (HSE06) for Hexi 
and -1.4eV (LDA)/-1.6eV (HSE06) for Heo. 

The He interstitial formation energy is only weakly 
affected by the treatment of exchange and correlation 



(LDA: 1.55eV, HSE06: 1.98eV). The same applies to 
binding energies. As in the case of alumina He inter- 
stitials are thermodynamically the most stable form of 
He. 



C. Yttria 

The defect formation energies for yttria shown in Fig. 3 
confirm the trends that were already observed for alu- 
mina and titania. The oxygen vacancy in yttria resem- 
bles its counterpart in alumina in so far as it also exhibits 
the characteristics of a deep defect. Helium interstitials 
occupy positions that are equivalent to 16c Wyckoff sites 
in the perfect lattice, which correspond to the structural 
vacancies of the bixbyite structure (see Sect. 3), and they 
bind to vacancies. 
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FIG. 4. Formation energies of He interstitials in various ox- 
ides as a function of the free volume at the interstitial site. 
Filled symbols indicate values at the respective equilibrium 
volumina. 



IV. HE INTERSTITIALS IN BULK OXIDES 

In the previous three sections it was demonstrated that 
within the thermodynamically allowed range of Fermi 
levels He interstitials prevail. The actual He intersti- 
tial formation energies for different oxides, however, vary 
over a wide range (AI2O3: 2.15eV, TiOii 1.55eV, Y2O3: 
0.71 eV, all values from conventional XC functionals). To 
resolve this variation He interstitial formation energies 
were computed as a function of volume for the oxides 
discussed above as well as the yttrium aluminum oxides 
and rocksalt structured oxides described in Sects. 4 and 5 
of the appendix. Since the XC functional was shown in 
Sect. HI to have a minor influence on He interstitial for- 
mation energies, calculations were carried out using con- 
ventional XC functionals that are computationally much 
more efficient than their hybrid relatives. 



A. Scaling relation for formation energies 

As shown in Fig. 4 the formation energies are found to 
scale remarkably well with the free volume at the inter- 
stitial site, where the latter is measured by the Voronoi 
volume of the He site in the relaxed configuration.'^^ It 
turns out that other measures such as the distance from 
the He site to the nearest neighbor atom (which is equiv- 
alent to constructing the smallest sphere around the He 
interstitial) do not yield such favorable scaling relations. 
Figure 4 also demonstrates that the formation energies 
of He interstitials in iron do not fall in the range of the 
oxides and for the same free volumes are systematically 
higher. 

The scaling relation in Fig. 4 is almost independent of 
the chemistry of the oxide involved and provides ample 
evidence that the major source of variation in the forma- 
tion energies of He interstitials in oxides is the volumetric 
compression of the He atom. It thus effectively decou- 



TABLE II. Migration barriers in eV for He interstitials in 
several oxides. Note that jump directions are approximate. 
AVrei' change in Voronoi volume of He site between initial 
state and saddle point normalized by volume of initial state; 
Avrei'- change in He-nearest neighbor distance between initial 
state and saddle point normalized by initial neighbor distance. 

Material Direction Barrier (eV) AVrei (%) Ar^ei (%) 





111 


2.16 


-12.5 


-11 




001 


3.86 


—9.2 


— 16 


TiOa 


001 


0.11 


3.6 


2 




110 


1.22 


-9.7 


-8 


Y2O3 


111 


0.70 


1.4 


-15 




111 


0.29 


4.0 


-9 




021 


2.73 


0.0 


-23 


YAP 


120 


1.22 


24.6 


-5 




101 


0.98 


3.6 


-3 


MgO 


100 


0.76 


16.1 


-10 


CaO 


100 


0.84 


17.2 


-9 


SrO 


100 


0.78 


17.4 


-9 


BaO 


100 


0.61 


17.6 


-8 
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FIG. 5. Migration barriers for He interstitials in several oxides 
as a function of the relative change in He-nearest neighbor 
separation. The inset shows the same data as a function of 
the relative change in the Voronoi volume of the He site. The 
light gray stripe is intended as a guide for the eye. 



pies chemistry from geometry and provides the basis for 
a simplified treatment of interfaces in Sect. V. 



B. Migration barriers 

To compare the mobility of He interstitials in oxides 
with iron their migration barriers in several oxides were 
determined. For each oxide we included paths that in- 
cluded up to three neighbor shells on the lattice of Wyck- 
off sites that correspond to the most stable He position. 
The results are compiled in Table II. 

To establish a scaling relation for migration barriers 
similar to the one for formation energies, we considered 
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FIG. 6. Binding energy of He interstitial clusters in several 
oxides as well as iron (Fe data from Ref. 5) as a function of (a) 
the number of He atoms in the clusters and (b) the density 
of He interstitial sites in the host material. 



the change in free volume (again measured using the 
Voronoi construction) from the initial He position Vini 
to the saddle point Vsad normalized with respect to the 
initial volume 



Vsad — Vi, 



(5) 



as well as the relative change in the distance between He 
and its nearest neighbor, Ar^e;, defined analogously. 

Figure 5 shows the migration barriers to scale fairly 
well with both of these measures, especially the relative 
He-nearest neighbor distance, but neither measure yields 
as convincing a scaling relation as in the case of the for- 
mation energies. The figure also contains the migration 
barrier for Fe,^ which at 60meV is substantially lower 
than any of the migration barriers found in the oxides. 
One can thus expect He interstitials in oxides to be much 
less mobile than in iron. 



C. Helium interstitial clustering 

It has been established by previous first-principles 
calculations^ that He interstitials in Fe exhibit a strong 
tendency to bind both to other He interstitials and vacan- 
cies. The resulting defect complexes represent nascent He 
bubbles and therefore play a key role in understanding 
He embrittlement in iron and steels. 

Complementing this information with regard to ODS 
steels. Fig. 6 shows binding energies of small He intersti- 
tial clusters in oxides. The left hand panel of Fig. 6 shows 
the dependence of the binding energy on the number of 
He atoms in the cluster for several oxides as well as iron. 
According to these data He-He interactions in oxides are 
typically much weaker than in iron with the exception of 
MgO. In some cases the interaction is even repulsive. 

The trends displayed in Fig. 6(a) can be readily ex- 
plained, at least qualitatively, by considering the density 



of He interstitial sites in a given host and thus the av- 
erage distance between nearest neighbor He interstitials. 
Since all He interstitials impose a relatively short-ranged 
compressive strain on the surrounding lattice, one can 
expect two He interstitials to interact weakly repulsively 
at moderate separations. In contrast, at very short dis- 
tances two He interstitials can effectively attract each 
other since the total lattice strain is confined to a smaller 
volume and thus the total strain energy is less for a He 
pair than for two separate interstitials. 

As shown in Fig. 6(b) the binding energies for small 
He clusters do indeed exhibit a pronounced dependence 
on the density of interstitial sites with larger densities 
being associated with stronger binding.'^* For iron, which 
exhibits a strong tendency to form He clusters but is not 
included in Fig. 6(b), the interstitial site density is 0.53 
per A"^ and thus much higher than in any of the oxides. '^^ 

The reasoning above and the data in Fig. 6 suggest that 
in general one can expect He clustering to be stronger in 
oxides with a higher He interstitial site density and that 
most oxides are less prone to He cluster formation than 
iron. 



V. OXIDE-IRON INTERFACES 

In the previous section properties of He defects in ox- 
ides have been extensively characterized, which led to the 
conclusion that He predominantly occurs in the form of 
interstitials. Furthermore it was found that formation 
energies of these defects scale with the free volume at 
the interstitial site almost independent of the chemistry 
of the host oxide. Similar though less pronounced trends 
could also be demonstrated for migration barriers and He 
interstitial cluster binding energies. 

The oxide particles in ODS steels vary widely in com- 
position and structure resulting in a rich variety of oxide- 
matrix interfaces of seemingly arbitrary complexity. Al- 
though some orientation relationships have been exper- 
imentally established (e.g., FejYAM, see Ref. 15), mod- 
eling these interfaces remains a computationally almost 
intractable task. The results of the previous section, how- 
ever, suggest that it should be possible to rationalize the 
behavior of He at oxide-iron interfaces primarily based 
on the geometry of the interface and the associated free 
volume. Based on this rationale we decided to focus on 
two model systems, FejMgO and Fe|FeO|MgO. These in- 
terfaces are particular suitable for a computational study 
since the lattice mismatch between Fe and MgO is small, 
which means that for the present purpose interface dis- 
locations can be neglected. 



A. Ideal interfaces 

The interface between two materials with rocksalt and 
body-centered cubic structure, respectively, is described 
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FIG. 7. (a,b) Illustration of the Baker-Nutting orientation 
relationship. Projection of (a) rocksalt MgO and (b) body- 
centered cubic Fe along [001] . The latter is rotated by 45 deg 
about the [001] axis such that the [lOOjrocksait and [110] Dec 
directions are parallel to each other, (c) Variation of interface 
energy as a function of the lateral displacement of the two 
crystals with respect to each other. 
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FIG. 8. Plane-averaged charge density and out-of-plane strain 
for the (a,b) Fe|MgO and (c,d) Fe|FeO|MgO interfaces as a 
function of position perpendicular to the interface plane. The 
colored spheres in the bottom panel indicate the atomic po- 
sitions. 



by the Baker-Nutting orientation relationship.'^^ It is il- 
lustrated in Fig. 7, which shows that the (OOl)rocksait and 
(OOl)bcc planes as well as the [lOOjrocksait and [110]bcc di- 
rections are parallel to each other. Among the earth al- 
kaline oxides MgO has the smallest lattice mismatch with 
Fe (4% based on the experimental lattice constants) and 
was therefore selected for the present study. Incidentally 
Fe-MgO interfaces have recently attracted a lot of at- 
tention since they exhibit a strong transverse magnetic 
resistance effect. ''-'^"''^ 

In the present work both Fe|MgO and Fe|FeO|MgO in- 
terface models were included, the geometries of which are 
illustrated in Fig. 8. First the minimum energy configura- 
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FIG. 9. Helium interstitial formation energies and their 
respective Voronoi volumina for the (a) Fe|MgO and (b) 
Fe|FeO|MgO interfaces as a function of position perpendicular 
to the interface. Helium atoms placed in the range indicated 
by the horizontal bars relax into the interface. 



tion of the Fe|MgO interface was established by scanning 
the energy landscape as a function of in-plane displace- 
ment, which yields the Fe-on-top-of-0 configuration as 
the most stable one. The Fe|FeO|MgO model is then ob- 
tained from the FejMgO model by inserting O atoms in 
the outermost Fe layers such that the Fe and O atoms 
form a square lattice parallel to the interface. Both mod- 
els were subsequently fully relaxed allowing both ionic 
motion as well as cell shape and volume changes. 

After relaxation of the Fe|MgO interface model the 
Fe slab is under a compressive in-plane strain of 4.8% 
while the MgO half is under a tensile in-plane strain of 
— 1.2%. The strain along [001] is of the respective op- 
posite sign near the interface as shown in Fig. 8(b) and 
quickly decays to the bulk value with increasing distance 
to the interface. For the Fe|FeO|MgO geometry the in- 
plane strains for the Fe and MgO part are 5.2% and ap- 
proximately zero, respectively. Distinctly non-zero strain 
along [001] is only observed for the Fe slab, for which the 
strain is tensile directly next to the FeO layer but com- 
pressive anywhere else. The different magnitudes of the 
strain in the Fe and MgO parts reflect the fact that the 
tetragonal shear modulus is almost 1.5 times higher for 
MgO than for Fe. 



B. Helium at interfaces 

Starting from the fully relaxed interface models. He in- 
terstitials were inserted sampling all distinct known bulk 
sites as well as sites in the interface. The calculated 
formation energies for these configurations are shown in 
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Fig. 9. 

For both interface configurations tlie formation energy 
for He interstitials is by far the lowest if the latter are 
located at the interface. For example in the case of the 
FcjMgO interface the formation energy is 2.7 eV to be 
compared with values of 3.3eV and 3.9 eV near the cen- 
ter of the MgO and Fe slabs, respectively. In the MgO 
part the formation energy is thus almost identical to the 
bulk value at the equilibrium lattice constant, whereas 
the corresponding value for Fe is noticeably lower than 
its unstrained bulk counterpart. This behavior is caused 
by two effects: Firstly, the tetragonal shear modulus is 
softer for iron than for the oxide, as a result of which the 
Fe slab is more severely strained. Secondly, the Poisson 
ratio of iron is less than 0.29, which implies that the av- 
erage volume per atom changes with strain. In the case 
at hand, one; observes an increase in the free volume at 
the He interstitial site from 7.4 A'^ to about 8.0 A^, which 
in combination with Fig. 4 readily explains the observed 
decrease in the formation energy.^^ 

Comparison of formation energies and Voronoi volu- 
mina shows that inside both the Fe and MgO slabs the 
formation energy scales with the free volume. At the 
FejMgO interface the lowest formation energy also cor- 
responds to the largest free volume. In the case of the 
Fe|FeO|MgO interface the formation energy is lower at 
the FeOjMgO interface than at the FeOjFe interface even 
though the free volumes suggest the opposite trend. This 
effect is related to the earlier observation that Fe does not 
fall on the same scaling relation as the oxides (see Fig. 4). 
While this implies that one cannot predict formation en- 
ergies at interfaces based exclusively on free volume, it 
nonetheless demonstrates clearly that He interstitials are 
more strongly attracted to interfaces than to either one 
of the bulk phases. For the two systems considered here 
as well as e.g., Fe|Y203 interfaces,^'' one observes an in- 
crease in the free volume at the interface compared to the 
bulk phases. This behavior can be attributed to a com- 
parably weak adhesion between the oxide and the iron 
matrix, which is related to the transition from mixed 
ionic-covalent to metallic bonding across the interface. 



VI. DISCUSSION AND CONCLUSIONS 

In this section the implications of the foregoing investi- 
gation for the understanding of He sequestration in ODS 
steels will be discussed. To simplify the discussion, first 
the major results of this work will be recapped. 

A detailed examination of three different oxides of dif- 
ferent chemical composition, stoichiometry and structure 
using both conventional as well as hybrid XC function- 
al showed that intrinsic defects limit the range within 
which the electron chemical potential can vary. Within 
the thermodynamically allowed range He interstitials are 
the most stable form of helium, yet they will bind to exist- 
ing vacancies. While conventional and hybrid functionals 
yield different values for band gaps as well as absolute de- 



fect formation energies for vacancies and substitutional 
He, the results are qualitatively similar. Furthermore, 
for He interstitials conventional and hybrid functionals 
provide formation energies that are also quantitatively 
similar. Since the latter are computationally much more 
demanding, further calculations employed conventional 
functionals only. 

Compared to iron He interstitial formation energies in 
oxides are significantly lower and exhibit a considerable 
materials dependence. Using data for a wide range of 
different oxides and volumes it was shown that the latter 
dependence can be described with good accuracy by a 
scaling relation based on the Voronoi volume of the He 
site. This finding demonstrates that the behavior of He 
can be largely rationalized in terms of free volume. It 
thereby greatly simplifies the task of understanding the 
behavior of Ho in ODS steels since it allows us to separate 
structure and chemistry. 

He interstitial migration barriers were found to be sys- 
tematically higher than in Fe. While the migration en- 
ergies do not obey a scaling relation as cleanly as for 
the formation energies, as a general trend the migration 
barriers decrease if the relative change in the He nearest 
neighbor distance from the initial configuration to the 
saddle point increases. 

In iron He interstitials bind strongly to each other, 
which facilitates the formation of He bubble nuclei.^ In 
contrast in most oxides the binding energies between He 
interstitials are much smaller indicating a weaker propen- 
sity for bubble formation. Binding between He intersti- 
tials was observed to scale with the volume density of 
He interstitial sites, with higher densities (e.g., MgO, Fe) 
leading to stronger binding. 

The variability and complexity of oxide-iron interfaces 
in ODS steels prevents a direct first-principles study of 
He sequestration at these interfaces. In the present study 
it was therefore decided to study two particular simple in- 
terfaces, Fe|MgO and Fe|FeO|MgO. In view of the scaling 
relations described above, one can expect these interfaces 
to act as prototypes for the types of interfaces that occur 
in real materials. The calculations revealed that in both 
interfaces the formation energies of He interstitials are 
significantly lower than in cither of the bulk phases. The 
low formation energies could again be correlated with a 
larger free volume at the interfaces. While in the MgO 
slabs the He interstitial formation energies reached al- 
most exactly the value that was obtained earlier for bulk 
MgO, formation energies for positions inside the Fe slabs 
deviated significantly from the bulk value. This behav- 
ior is directly related to strain fields that distort the Fe 
lattice and thereby affect the free volume available at in- 
terstitial sites. The strain effect is much larger in the Fe 
matrix because its tetragonal shear modulus is consider- 
ably lower than the one of MgO. 

Now one can combine all this information to develop 
a schematic energy landscape for He interstitials in ODS 
steels. In the Fe matrix He interstitial formation ener- 
gies are high but migration barriers are low. In contrast. 
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FIG. 10. Schematic energy landscape for He interstitial mi- 
gration in an ODS steel. In the Fe matrix formation energies 
are high but migration barriers are low, while the opposite 
applies for the oxide particles. The smallest formation en- 
ergies and thus the highest solubilities are predicted in the 
interface region. Strain fields can lead to gradients near the 
interface that depending on the sign of the strain field can 
either increase or decrease toward the interface. 



formation energies in the oxide particles are lower while 
migration barriers are higher than in Fe. The lowest 
formation energies are observed at the interface. Since 
iron is elastically softer than the majority of oxides con- 
sidered in this study, the iron matrix is more likely to 
be strained. This affects the free volume available to 
He interstitials and accordingly their formation energies 
(and solubilities). Combining these data, one can ob- 
tain a schematic energy landscape as the one sketched 
in Fig. 10. In this particular plot the free volume in 
the Fe matrix is assumed to increase toward the inter- 
face leading to a gradual sloping of the landscape toward 
the oxide particle. It is, however, equally possible that 
the free volume decreases toward the interface, leading 
to a landscape that rises toward the interface. In real- 
ity the character of the strain field sensitively depends 
on the structure of the interface, which can include e.g., 
interface dislocations, amorphous regions, and chemical 
gradients (in so far as they translate to strain). In fact, 
experimentally He bubble formation is observed not to 
occur around all precipitates,^^ for which strain effects 
could be a reasonable explanation. 

In the introduction, it was pointed out that oxide par- 
ticles in ODS steels exhibit a broad spectrum of struc- 
tural and chemical variations, including amorphous re- 
gions and extremely small nanometer-sized inclusions. 
While these structures were not explicitly studied here, 
the present work has established clear trends that hold for 
a variety of different local environments both structurally 
and chemically. It was also demonstrated that He in- 
corporation in oxides is qualitatively different from iron. 
This is most clearly visible in the He interstitial formation 
energies (Fig. 4), for which oxides display systematically 
lower values than iron for the same free volumes. This 
behavior can be readily understood in terms of mixed 
covalent-ionic bonding in the oxides vs metallic bonding 
in iron. The covalent character of oxides leads to stronger 



directional bonding and a more localized electron charge 
density compared to iron. As a result oxides adopt more 
open structures with larger interstitial sites (both in ionic 
and electronic terms), and defect induced strain fields are 
less extended than in metals. These qualitative features 
are also found in amorphous oxide particles and small 
oxide inclusions. (As shown in Refs. 10 and 12, already 
very small oxide clusters containing just a few atoms fea- 
ture pronounced directional bonding). It is therefore ex- 
pected that the results obtained in the present work also 
transpire to these more general situations. 

In summary, in this paper it has been demonstrated 
that oxide particles in ODS steels have a higher He solu- 
bility than the Fe matrix, which is primarily the result of 
larger interstitial sites. The solubility at oxide-iron inter- 
faces is even larger than in the bulk oxides. Strain fields, 
which can affect in particular the iron matrix, lead to 
solubility gradients near oxide-iron interfaces. The data 
obtained in this study not only provides valuable insight 
into the behavior of He in ODS steels but can be used to 
derive parameters for rate equation models of He seques- 
tration in ODS steels.^'' 
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Appendix: Bulk properties of oxides 

This appendix summarizes results for groundstate 
properties of the oxides included in this study as obtained 
using different computational methods. Based on these 
results the parameters for the defect calculations were 
chosen, in particular the XC functionals. 

1. Alumina 

The ground state structure of alumina (AI2O3) is 
corundum (Strukturbericht symbol D5i , space group R3c, 
number 167), which has a primitive unit cell of rhombo- 
hedral symmetry containing ten atoms with Al and O 
atoms occupying Wyckoff sites 4c and 6e, respectively. 
The structure can also be described using a hexagonal 
setting, in which case the unit cell compromises three 
times as many atoms. The relation between the two set- 
tings is discussed in detail in Ref. 55. Following com- 
mon procedure the structural parameters in Table HI 
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TABLE III. Structural and electronic properties of AI2O3, 
Ti02, and Y2O3 from experiment and calculations, a, c: lat- 
tice parameters (A), 9: rhombohedral angle (dcg), Xi, yi, Zi, u: 
internal structural parameters, Eg'- band gap (cV), epsilorioo- 
electronic contribution to dielectric constant, epsiloriion- ionic 
contribution to dielectric constant, epsilono: static dielectric 
constant (sum of £00 and eion). The subscripts -L and || indi- 
cate the dielectric constant perpendicular and parallel to the 
rlionibolicdral [111] axis (equivalent to the [0001] axis in the 
hexagonal setting), respectively. Note that for aluniiiia the 
thermal band gap is reduced with respect to the optical gap 
(given here) due to polaronic effects (Ref. 48). Experimen- 
tal data for alumina from Refs. 48 and 49, for titania from 
Refs. 50-52, and for yttria from Refs. 53 and 54. 

Alumina (corundum, AI2O3), space group R3c, number 167 
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0.3523 


0.3520 


0.3522 


0.3519 
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Yttria (bixbyite, Y2O3), space group lo3, number 206 




Expt. 
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are given in the hexagonal setting. The data shows that 
LDA calculations provide slightly better agreement with 
structural parameters than PBE calculations. The band 
gap underestimation in these calculations is typical for 
conventional XC functionals. This deficiency is largely 
corrected by using hybrid functionals, which also further 



improve the structural parameters. For defect calcula- 
tions the PBE and HSE06 functionals were used. 



2. Titania 

There is a wide range of titaniiim oxides with differ- 
ent stoichiometrics and crystal structures. Here the ru- 
tile structure {Strukturhericht symbol C4, space group 
P42/mnm, number 136) is considered, which has a prim- 
itive unit cell of tetragonal symmetry containing 2 Ti 
and 4 O atoms, which occupy Wyckoff sites 2a and 4/, 
respectively. For the structural parameters the closest 
agreement with experiment is obtained for the LDA and 
HSE06 functionals, which were accordingly used in defect 
calculations. 

The PBE functional yields extremely soft phonon 
modes at F, which leads to a mimerical divergence in 
the ionic dielectric constant. In reality Ti02 does indeed 
possess soft modes that give rise to a very large dielectric 
constant. Due to their softness it is not trivial to accu- 
rately reproduce their frequencies and more crucially the 
resulting ionic contribution to the dielectric constant. In 
fact, the PBE functional did not only caused problems in 
calculations of the dielectric response but also in a num- 
ber of defect calculations involving He interstitials: After 
removing the He interstitial from the relaxed supercell, 
the configuration did not relax to the original lattice po- 
sitions but instead assumed a distorted striicture with an 
energy that was lower than for the perfect rutile lattice. 
These problems raise concerns regarding the reliability 
of the PBE functional for calculations in Ti02. These 
problems were specific for the PBE functional and did 
not transpire to the HSE06 hybrid functional. 

3. Yttria 

The ground state structure of yttria (Y2O3) is bixbyite 
{Strukturhericht symbol DSs, space group Ia3, number 
206). The primitive cell contains 16 Y atoms on Wyck- 
off sites 86 and 24rf as well as 24 O atoms on Wyckoff 
sites 48e. The bixbyite structure can be described as a 
calcium fluorite lattice in which one quarter of the an- 
ion sites are unoccupied ("structural vacancies"). Struc- 
tural and selected electronic properties obtained using 
different XC approximations are shown in Table III in 
comparison with experimental data demonstrating good 
overall agreement. For this material defect calculations 
were carried out using the PBE XC functional only. 

4. Yttrium aluminum oxides 

There are three distinct ground state structures in the 
Y-Al-0 system that contain all three elements, a per- 

ovskite (yttrium aluminate, YAP), a garnet (YAG), and 
a monoclinic (YAM) structure. These ternary oxides are 
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TABLE IV. Structural and electronic properties of yttrium 
aluminum oxides. Symbols as in Table III. Experimental 
data for YAP, YAG and YAM from Refs. 56, 57, and 58, 
respectively. 

YAP (YAIO3), space group Prima, number 62 



Experiment LDA PBE 



a 5.315 


5.284 


5.392 


b 7.354 


7.314 


7.452 


c 5.167 


5.133 


5.226 


XYi 0.053 


0.054 


0.055 


ZYi 0.988 


0.987 


0.987 


xoi 0.477 


0.480 


0.477 


zoi 0.083 


0.083 


0.086 


2:02 0.292 


0.293 


0.294 


yo2 0.042 


0.044 


0.045 


Z02 0.707 


0.706 


0.706 


YAG (Y3AI5O12 


), space group Io3d, number 230 


Experiment 


LDA 
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12.006 


11.916 
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xo -0.032 


-0.031 


-0.031 


yo 0.051 


0.051 
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0.149 
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YAM (Y4AI2O9), space group P2i/c, number 14 


Experiment 


LDA 


PBE 


a 7.375 


7.296 


7.441 


b 10.462 


10.380 


10.563 


c 11.110 


11.044 


11.260 


/3 108.58 


108.46 


108.50 



included because yttrium and aluminum oxide are com- 
mon additions in ODS steels and some of the resulting 
oxides have been observed experimentally.^^ 

The perovskite structure (YAP) has orthorhombic 
symmetry with space group 62 (Prima). Its primitive 
unit cell contains 4 Al, 4 Y, and 12 O atoms correspond- 
ing to a chemical sum formula of YAIO3. Yttrium occu- 
pies 4c Wyckoff sites with two internal parameters, Xy 
and zy while the Al atoms sit on the 46 Wyckoff sites 
with no internal parameters. The oxygen atoms occupy 
both 4c and 8d sites. Their positions can be described 
using five internal parameters, xqi, zqi, X02, yo2, and 
zo2- Table IV compares the structural data obtained 
from experiments and calculations. 

The garnet structure (YAG) has cubic symmetry with 
space group 230 (Ia3d). Its primitive unit cell contains 
20 Al, 12 Y, and 48 O atoms corresponding to a chemi- 
cal sum formula of Y3Al50i2. Yttrium occupies the 24c 
Wyckoff sites and aluminum atoms are located both on 
16a and 24(i Wyckoff sites. Oxygen atoms occupy the 96h 
Wyckoff sites with three internal parameters, xo,yo, and 
Zo- The structural data obtained from experiments and 
calculations are compared in Table IV. 

Monoclinic yttrium aluminum oxide (YAM) belongs to 
space group 14 (P2i/c). Its primitive unit cell contains 
8 Al, 16 Y, and 36 O atoms corresponding to a chemi- 
cal sum formula of Y4AI2O9. Yttrium, aluminum, and 



oxygen occupy 4, 2, and 9 different 4e Wyckoff sites, re- 
spectively, leading to a total of 45 internal parameters. 
Table IV compares the lattice constants to experimental 
data revealing good overall agreement with the reference 
data. We refrain from listing the internal structural pa- 
rameters but assert that they are also in good agreement 
with experiment. 



5. Oxides in the rocksalt structure 

There are a number of oxides that adopt the rock- 
salt structure. Here only the earth alkaline oxides are 
included as the most simple variant. One obtains lat- 
tice constants of 4.240 A (experimental value 4.207 A, 
Ref. 60), 4.832A (4.803A, Ref. 60), 5.207A (5.160A, 
Ref. 61), and 5.611 A (5.524 A, Ref. 62) for MgO, CaO, 
SrO, and BaO, respectively. 
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